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© Verfahren und Schaltungsanordnung zum Ermitteln des Reibwerts. 



© Zum Ermitteln des Reibwerts zwischen einem Kraftfahrzeugreifen und der Oberflache einer Fahrbahn mit 
einer Schaltungsanordnung (1) werden die Fahrzeuggeschwindigkeit, die Drehgeschwindigkeit eines Rades (2-5) 
und der Radbremsdruck (p) laufend gemessen. Aus der Drehgeschwindigkeit wird die Winkelbeschleunigung 
des Rades berechnet und aus der Winkelbeschleunigung und dem Bremsdruck wird der Reibwert (u) mit einem 
rekursiven Schatzalgorithmus radselektiv ermittelt. ZweckmaBig 1st es, einen Schatzalgorithmus nach der 
Methode der kleinsten Quadrate zu verwenden. Der als Funktion des Radschlupfes ermittelte Reibwert dient zum 
genaueren Regeln eines Antiblockiersystems (10) oder eines Schlupfregelsystems. 
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0, Erfindung be* * EMMS 'SSSS^Sf 

ch. 1 bzw. 7. Sie d,en, znm ^™ I? R e jb w grts kdnnen Antiblockiersysteme besonders 

KSl to „ di. genau. Ke„n t „,s da. ™^;°„~ n den ^k* das Leak— 
Be, .,nem bekann,en R*^ ^^Tl^^^L^ und lies, aes .in., 
und die Fahrzaaggeschwindigkeit. Bn Pmzessor vmroe Lenkmoment zu 

zztjsz — CTdTB^:rr «». - ^ . *- 

mi, eine, Fermel berechnet, di. .en eine m K, al .eo.e,ehg.w,eh, »«lTi der Hinterrader 
di . Signal, von Sense,.. fO, die ^^TTi^^^^^ * f^™* 

Reibwert fur eine Antiblockierregelung zu ve .7 e " d n en vergnder|ichen Beibwert zwisc hen dem Kraftfahrzeugrei- 
Der Erfindung liegt die Aufgabe zugrunde, den ve ^rlichen ™ D s0 „ eine 

, fen und der Fahrbahn mit geringem Aufwand moghchst genau zu bertmmen. 
verbesserte Radbrems- und Radschlupfregelung ermoglicnen. . nn . 1 und 

die Schaltungsanord- 

niedergelegt. . _ H Rp jhwfirt mit bewahrten und effektiven 

' ^^TZtr.TTZT^r^'Z de, —and end der 
Speicherbedarf in Grenzen halten. . . . _ anhanr , ripr zeichnunq erlautert. Es zeigen: 

Fiaur 3 ?n m eSE!SModell fUr die erfindungsgemaSe Reibwerterrnittiung^ 

Abhangigked von ^^^r^fL, „ „ eist 2 „.i Vord.rrad.r 2 end 3 sowie zwei 
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di. Ob. S,.n Bus 32 d.m V.ndlb.eck 10 zugerobr, ™,den Die.., Ohrt. da™*,n £ e^-n^R* ^ 
und seine Lenkfahigkeit erhalten bleibt. 
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(1) 



dem Material, der Struktur und dem Zustand de Fa^Z 2ST p ^ im wesentlic ^n von 
z.B. zwischen einer trockenen Asphaltfahribahn und ein« n* variiert in ? inem weiten Bereich, 

hier auch radselektiv, d.h. fUr jedes hS^T^Z^T^ **■ *~ & wird 

zeugs bekanntlich auf Fahrbahnstellen n^T^^^'S^^ ^ 6ineS Kraftfahr ' 
anderen auch die Laufflachen der einzelnen Relfen^en^™ h- ^' 9to,t befmden k6nnen und da ™ 
Eine Reibkraft kann nur durch einen RadscWupfT UnterSCh,9dl,chen 2usta "d aufweisen konnen. 



25 Punkt u = 0 fur s = 0 unS steigen zu ■ dar9este,lten Kufven beoinnen bei dem 

von etwa s = 0.005...0 2 liegT Danach nimm £,?r ^ T*™™ m " daS bei einem ^schlupf 

KraftUbertragung wesen.chentaf^ ^ * Ub ^ 9 ™ - Or die 

« H 6 A d p;? t dVerZ09erU ; 9 6rkannt W6rden - Bn StarkeS Anstei 9 en der Ra^rzflflemniSJrSS^S 
ih^ t h r rQerat Bremsdruck P verrin 9ert. Wenn daraufhin cfc-Radw^B^XlTZS^Jf 
erhoht das Steuergerat wieder den Bremsdruck. Die sich dabei ergebenden zyklischen 1^11^ 

ISZXeT. V h'' 5 p^"^" An «blockier Sys teme 9 sind nach En^n£ZS 

ausgelegt. Die h.er beschnebene Ermittlung der Reibcharakteristik u = w(s) , d .h. der Abhangigkeit des 
Re.bwertes von dem Radschiupf, ermoglicht es, da sie in Echtzeit durchgefUhrt werden kann i 
Antiblockierregelung zu systematisieren. ' e 

45 er S ic D hLh r S ttl 7 Un9 t 0de h Sch H atZU " 9 des Reibwerts basiert a "f ^nem Ein-Rad-Modell, das aus Figur 3 
ers.chtfich .St. Zusatzhch zu den bereits erwahnten GroBen sind in dieser Figur noch dargestellt: der Anteil 

p a A J er / p ZeU9ma ' Se ' V °" dieS6m ein9n Rad 9 etra 9 en wird - das Massentragheitsmoment 9„ des 
Rades, das Bremsradius r B und die Bremskraft F B . 

Die Reibkraft F R verzogert die Fahrzeuggeschwindigkeit wie folgt: 

~ m A*A = F R , (3) 

55 wobei die Differenz zwischen dem Brems- und dem Reibungsmoment die Raddrehung wie folgt verzogert: 
G#>r = F B .r B - F R >r R (4) 
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Die Reibkraft ergibt sich aus: 
F R = u(s)F z = u(s)'m A -g (5) 
5 Mit einem Reibwert U B an der Bremsscheibe oder -backe kann die Bremskraft wie folgt ausgedrUckt 



werden: 



F b = Ub'Pb'A b (6) 

,o worin p s der Bremsdruck und A B die wirksame Druckflache des Bremsbelages aufder Bremssscheibe oder 
-trommel ist. Somit ergeben sich folgende Systemgleichungen fUr das Ein-Rad-Modell: 

Die Radgeschwindigkeiten werden Ublicherweise durch aufeinanderfolgende induktive Abtastwerte gemes- 
20 sen. Die Zeitinterval.e zwischen den Abtastwerten verandern sich mit der Drehge schwmd.gk a t Aus 
Vereinfachungsgrunden wird hier eine konstante Abtastzeit T A = 10 ms gewahlt. Dam.t erg.bt s.ch folgende 
diskrete Beschreibung des Ein-Rad-Modelis. 

co , (") - a> j. (w - 1) , VjMlL p (») - ^gg- fi(n) (8) 
T A & R " ® R 

Die Raddrehgeschwindigkeit Mn) und der Bremsdruck p B (A» werden I laufcnd ge™ssen ^ " der 

30 Zahlindex fur die aufeinanderfolgenden Messungen. Der tatsachl.che Re.bwert u(n) wd njejem 
rekursiven Schatzalgorithmus nach der Methode der kleinsten Quadrate - in der L.teratur als RLS-Verfahren 
beschrieben - ermittelt. Mit einer MeBvariable y(n) 

y{n) = ~ 0 + ^£L pB jn) (9) 

T A ©„ 

und einem Fahrzeugparameter X 

40 

x = Ta3Ie. {10 ) 

45 

ergibt sich folgender rekursiver Schatzalgorithmus: 
u(n) = E(n-1 ) + /f(n)(y(n)-ru(f7-l)) f/ /j 
50 mit einem Fehlerruckfuhrungs-Verstarkungsfaktor k(n) 

*w .,A-i) x ' f< "" 1) ' 

55 

worin die P(n-1) die Kovarianz des Schatzfehlers ist 
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P(/)-1) = 10^(y(A7-1)-A > u(n-1))2. (13) 




nW = q(i-«*')-c,j 



(14) 



15 




U(S) = AS + &JS, (15) 

die folgenden linearen Schatzalgorithmus ergibt: 

u(n)=As(n)+By/s(n). (16), 

25 worin A und B Parameter einer Kurve sind, mit welcher der zu ermittelte Reibwert approximiert wird Sie 
werden von den als Vektor 

b T =[AB] 

30 dargestellten Kurvenparametern A und B, von den als Datenvektor 



dargestellten MeBwerten des Radschlupfes s und von dem wie weiter oben beschrieben ermittelten 
35 Reibwert u(n) abgeleitet. Die Berechnung des Pradiktionswertes fur die Reibcharakteristik wird mit nachfol- 
genden Gleichungen durchgefuhrt 

u{n) = X r {n)'b(n) 



k{n) = P(n-1 )X(/i)(1 + X T (n)P(n--\ )X(n))~' 

45 Ein Schatzwert u(n) fur den Pradiktionszeitpunkt n wird aus den MeGergebnissen zu dem Zeitpunkt n-1 
ermittelt. Ein neuer Schatzwert 



zum Zeitpunkt n wird entsprechend aus dem Schatzwert zu dem Zeitpunkt n-1 berechnet. Ein Vektor k(n) 
zum Zeitpunkt n wird ebenfalls aus den Werten zum Zeitpunkt n-1 berechnet. Dieser Vektor wird in der 
Literatur auch als Kalmanscher Verstarkungsfaktor bezeichnet. Die in ihm enthaltene Kovarianz-Matrix P(n- 
1) wird mit der Gleichung 



X T =[sVs] 



40 



b{n) = b(n - 1) + £(«)(u.(k) - X T (n)b(n - 1)) 



bin) 



55 



P(n)= \0" \n{n)-X r (n)b(n)\ 2 -l 



(18) 
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berechnet, allerdings mit den Abtastwerten vom Zeitpunkt n-1. Sie hangt vom Schatzfehler ab. I ist hier die 
Einheitsmatrix. 

Das hier zugrunde gelegte Modell gemaB Gleichung (15) fUhrt - wie auch aus dem Diagramm von Figur 
5 ersichtlich ist - zu guten Ergebnissen, wenn Einschrankungen fUr den Ort des Kurvenmaximums 
zugelassen werden. Die groBeren Abweichungen nach dem Maximum konnen toleriert werden, da moderne 
Antiblockiersysteme in diesem Bereich nicht arbeiten. 

Ein wahlweise hier verwendetes Modell der Reibcharakteristik 

M = ~ — 7TT (19) 

verwendet drei Parameter A, B und C, die eine gute Approximation fur den Reibwert Uber den ganzen 
Radschlupfbereich ergeben. Die Ausgangssteigung der Kurve u(s) ist 

u(0)=1/>4 (20) 

und das Maximum des Reibwerts betragt 

(21) 



und zwar bei einem Wert s 0 des Radschlupfes 
S 0 =V^7C- (22) 

Ein Modell der Reibcharakteristik nach einer Beziehung 

y (s)=s=A+Bs + Cs 2 



ergibt ebenfalls einen linearen Schatzalgorithmus. 

Der beschriebene Schatzalgorithmus minimiert das Fehlerquadrat des Schatzfehlers. Vorteilhaft ist, dafl 
alle bisherigen MeBwerte in die aktuelle Schatzung eingehen, womit fur die Schatzung mehr Information zur 
Verfugung steht. Dadurch konnen Fehler der Radaufstandskraftschatzung besser als mit bekannten Verfah- 
ren kompensiert werden. 

Patentanspriiche ' 

1. Verfahren zum Ermitteln des Reibwerts zwischen einem Kraftfahrzeugreifen und der Oberflache einer 
Fahrbahn, bei dem die Fahrzeuggeschwindigkeit und ein weiterer Fahrparameter gemessen und aus 
diesen MeQwerten der Reibwert abgeleitet wird, 

dadurch gekennzeichnet, 

- daB die Drehgeschwindigkeit eines Rades und der Bremsdruck des Kraftfahrzeugs laufend 
gemessen werden, 

- daB aus der Drehgeschwindigkeit die Winkelbeschleunigung des Rades berechnet wird, 

- und daB aus der Winkelbeschleunigung und dem Bremsdruck der Reibwert (u) mit einem 
rekursiven Schatzalgorithmus ermittelt wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB der Reibwert (u) mit einem rekursiven 
Schatzalgorithmus nach der Methode der kleinsten Quadrate (RLS-Verfahren) radselektiv ermittelt wird. 

3. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB der Reibwert (u.) mit folgendem rekursiven 
Schatzalgorithmus ermittelt wird: 

U(n) = u(n-1 ) + /f(n)(y(/7)-X*il(n-1 )), worin 
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eine aus den Mefiwerten gebildete Variable, 



x= m * gr * 



ein Fahrzeugparameter, 



*(»).- 



l + X 2 . 



'P(n - 1) 



X * P(n - 1) 



ein Verstarkungsfaktor fUr die FehlerrUckfUhrung und 

P(n-1) = 10"*(y(n-1)-X*Il(n-1)) 2 die Kovarianz des Schatzfehlers zum MeBzeitpunkt n-1 sind. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB der Reibwert (u(s)) als eine Funktion des 
Radschlupfes (s) mit einem rekursiven Schatzalgorithmus ermittelt wird. 

Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB der Reibwert auf der Grundlage eines 
Modells U(s) = As + By/s ermittelt wird, worin A und B Kurvenparameter mit der Vektordarstellung 



und aus dem ermittelten Reibwert v(n - 1) berechnet werden, und daB ein neuer Reibwert u(n) mit 
folgendem Schatzalgorithmus ermittelt wird: 
n(n) = X T (n)'b(n) , worin 



b T = [AB] 



sind, die aus den gemessenen Radschlupfwerten mit der Vektordarstellung 



X T = [sJs] 



b{n) = b{n - 1) + k{n - X 1 \n)b(n - 1)) 



k(n) = P(n-1)X(n)(l + X T (n)P(n--\)X(n))-i und 



P(n - 1) = 1 0' * - 1) - X T (« - l)b(" ~ 1)| 2 



*/ sind. 



Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB der Reibwert (u(s)) auf der Grundlage eines 
Modells 



s 



H(5) = 



A + Bs + Cs 1 



ermittelt wird, 
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7. Verfahren nach Anspruch 4, dadurch gekennzeichnet, daB der ermittelte Reibwert (u) zum Steuern 
eines Antiblockiersystems (ABS) oder eine's Schlupfregelsystems verwendet wird. 

8. Schaltungsanordnung zum Ermitteln des Reibwerts zwischen einem Kraftfahrzeugreifen und der Ober- 
flache einer Fahrbahn gemaS einem Verfahren nach einem der vorstehenden Anspruche, die versehen 
ist mit Sensoren (24), durch welche die Fahrzeuggeschwindigkeit und ein weiterer Fahrparameter erfaBt 
werden, sowie mit einer Auswerteschaltung (20), durch welche die Sensorsignale ausgewertet und aus 
ihnen der Reibwert berechnet wird, dadurch gekennzeichnet, 

- daB durch die Sensoren (24) die Drehgeschwindigkeit eines Rades (2-5) und der Bremsdruck (p) 
an einem Rad laufend gemessen werden, und 

- daB durch die Auswerteschaltung (20) 

-- aus der Drehgeschwindigkeit die Winkelbeschleunigung des Rades (2-5) berechnet und 

~ aus der Winkelbeschleunigung und dem Bremsdruck der Reibwert (u) mit einem rekursiven 

Schatzalgorithmus ermittelt wird. 

9. Schaltungsanordnung nach Anspruch 7, dadurch gekennzeichnet, daB der Reibwert (u) mit einem 
rekursiven Schatzalgorithmus nach der Methode der kleinsten Quadrate (RLS-Verfahren) ermittelt wird. 

10. Schaltungsanordnung nach Anspruch 1, dadurch gekennzeichnet, daB sie zum Regeln eines Antiblok- 
kiersystems (ABS) Oder eines Schlupfregelsystems verwendet wird. 
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